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Abstract 

Introduction. The movement of the carrier and external factors (the effects of the atmosphere, temperature and 
pressure) degrade significantly the image quality of the servo optoelectronic systems (OES) and the positioning 
accuracy of the emitting OES. The issues of image quality improvement and the probability of keeping the image of the 
observation object (OO) on the optical axis of the servo EOS are considered. 

Materials and Methods. The development of an automatic control system for an optoelectronic device (ACS OED) 
involved solving a multi-criteria optimization problem taking into account a number of conflicting technical-and- 
economic (TE) requirements. The determination of tolerated dynamic errors (TDE) of image stabilization was a key 
issue in the development of on-board optoelectronic devices (OOED). Lagrange equations of the second kind and the 
mixed Gilbert method made it possible to obtain a mathematical model of the CO OED. Then, the decomposition of a 
two-link ACS with nonlinear cross-couplings in the CO was performed. A functional diagram of the image formation 
model of the OOED was presented. The parameters of the matrix photodetector and the requirements for the dynamic 
error of the ACS OED, taking into account the permissible MTF of the OED, were listed. The functions of transferring 
modulation, as well as linear, harmonic and vibrational shift of the image corresponding to the permissible and achieved 
TDE were visualized. Logarithmic frequency characteristics were created in the Mathcad environment. The two-link 
control system of the OED with the specified parameters of the CO for the considered movement was presented as two 
independent azimuth and elevation control channels. 

Results. The processes of control of the on-board optoelectronic system in the stabilization and tracking modes were 
described. To study the dynamics of spatial control of the OOEP in accordance with the ACS methodology, a computer 
simulation model (CSM) of the digital automatic control systems (DACS) of the OED was developed. It was 
implemented in the Matlab environment and consisted of CSM CO, drives, proportional-integral-derivative (PID) 
controllers taking into account non-linearities, a central computing device (CCD), a guidance software device, a CSM- 
carrier that implemented the equations of motion. Harmonic vibrations of the carrier were described. The errors of 
tracking and stabilization in the tracking mode with an additional control action introduced in the form of a constant 
speed were determined. The dynamics of spatial control of the OOED was investigated. A computer simulation model 
of a digital automatic control system of an optoelectronic device, the results of modeling the DACS OED without 
considering the board movement, and the processes of OED control subject to movement were visualized. 

Discussion and Conclusions. The stabilization accuracy was calculated for the studied cases. It was established that the 
stabilization tens of times exceeded the previously stated indicators, and it tens of times reduced the requirements for 


the convergence of the laser beam and the laser radiation power when developing the optical path of the product in 
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question. The proposed CSM can be used in the development of the on-board optoelectronic systems. In this case, the 
application of the presented methodology and CSM will help to reduce labor costs and minimize errors. 
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Introduction. When designing modern on-board optoelectronic devices (OED) and complexes, methods of 
computer modeling of optoelectronic systems (OES) are widely used. Such authors as Yu. G. Yakushenkov, 
V. V. Tarasov, I. P. Torshina, V. P. Ivanov, V. A. Baloev, V. A. Ovsyannikov, V. L. Filippov, worked on these tasks. 
Computer modeling enables to solve the problems of rational choice of the structure, parameters, and element base of 
the OES, providing the required performance indicators under specified constraints without expensive field studies and 
tests. 

The movement of the carrier and external factors (temperature, pressure) degrade significantly the image quality of 
the servo OES and the positioning accuracy of the emitting OES. The probability of holding the image of the object of 
observation (OO) on the optical axis of the tracking OES is reduced. For emitting OES, the positioning accuracy and the 
probability of performing observation tasks are reduced. To hold the image in vision or in the irradiation sector, it is 
required to capture a fairly large solid viewing angle. Hence, it is required to increase the dimensions of optical systems 
and power electronics capacities. To avoid this, controlled radiating OES is used. The above-mentioned features of the 
OED limit the development of the ACS OED. This problem was considered by V. A. Strezhnev, V. M. Matrosov, 
A. S. Zemlyakov, N.N.Malivanov, E.I.Somov, A.I. Malikov, V.A.Krenev, A.I. Karpov, D. A. Molin, 
A. V. Mikhalitsyn. When designing controlled OED operating in guidance and tracking modes, the following fact was 
established: the guidance time and the accuracy of stabilization of the optical axis, as well as the dynamics of the OED 
subsystems affect significantly the image quality. 

Thus, the task of minimizing the time and increasing the accuracy of pointing and holding the image of the OO in 
vision of the controlled OED is a challenge. In this regard, the following is required: 

— development of mathematical models of power lines as control objects; 

— synthesis of control algorithms; 

— creation of computer simulation models (CSM); 

— study on control systems that take into account the dynamics of the movement of a controlled OED; 

— determination of parameters affecting the dynamic properties and image quality. 

The work aimed at improving the accuracy characteristics and image quality of OED operating in the guidance, 
stabilization and tracking modes, due to the rational choice of their parameters during synthesis and modeling. 

Materials and Methods. The development of ACS OED started with solving a multi-criteria optimization problem 
that took into account conflicting technical-and-economic (TE) requirements. In the development of study [1], a 
modified method of designing a self-propelled ACS OED was proposed. For the quality criterion of the ACS, we took a 
set of dynamic characteristics of control channels that met the conditions: 

Toon (Vy) = Teen (Vy), 8(@) < 64°" (w), Aa, < Aa", Ad, < Ade", (1) 
M < M*°" |A@| = Ag?”, |AL| => AL*°". 

Here, T3" (Vy), Toon(Vu.) — permissible and implemented modulation-transfer function (MTF) of the OED at the 
Nyquist frequency, 5(@); 5°" — amplitude frequency response of cross-couplings and its acceptable decomposition value; 
Aol", Adio", Aa,, Ad, — permissible and steady-state values of the dynamic errors of the ACS in angle and angular 
velocity under the action of disturbances under conditions close to the actual operation of the ACS; k= 1,2, ... — 
number of the control channel providing image quality; Myo, = 1.05—1.25 — oscillation index, Ag; AL — stability 
margins in phase and amplitude. 

The key issue in the development of the OED is the determination of tolerated dynamic errors (TDE) of image 
stabilization. 


The image formation scheme of the on-board optoelectronic device (OOED) is shown in Figure 1. 
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or Carrier model 


ACS model 
ISS model 
. device 
OO model MPD model ACC model model 


Fig. 1. Functional diagram of the OOED imaging model: OO — object of observation; OS — optical system; MPD — matrix 
photodetector; ACD — amplifying-converting device; ISS — image stabilization system 


The MTF OOED should meet the condition that provides acceptable image quality [2]: 
Tosc (v) = Tay (V)To6 W)T on W)Tyny (v)Tccu (v) > To3c (v), (2) 


(mi) | 


-1 
om = . Yu 
To5c() = 2 (sinc (w is) exp - EGa 
Here, To3c(v) — MTF OED; N — spatial frequency; T35¢(v) — admissible MTF OOED; T,,,(v) — atmospheric MTF; 
Too(v) — lens MTF; Ty,,(v) —photodetector MTF; Ty,y(v) — optical information conversion function MTF [1]; 


Tccu(v) — MTF of the image shift (dynamic error of the ISS), depending on the type of dynamic displacement of the 
image: linear (JI) — x(t) = Vt, harmonic (() — x(t) = agsin(t), random (R), and defocusing. 
The permissible MTF of the image stabilization system includes tracking, stabilization, vibration protection (VPS) 


and automatic focusing (AFS) systems [3]: 


Ta N 
Tocu(N) = Tr (N)T-(N)Tg(N) Ta (N) = TEEN) _ oe 
Ty (v) = sinc(nAay,v),T-(v) = Jo(2tAdopv), Tg(v) = exp[—2(nAagv)?], as 
2),(A he ie 
To(v) = “Ate ay = 2810 (~*) (: _ =) 


Tor(N) = Tar(N )To6(N) Ton (NT yny CN). 
Here, T;(N) — image linear shift error; T-(V) — sinusoidal error amplitude; Tg,(N) — average value of vibration 
amplitude error; J, — zero-order Bessel function of the first kind; J, (Ay) — Bessel function of the first kind of the first 
order; Ay — focus error (mm); 6 — average value of wave aberration in fractions of a wavelength; 4,, — average 


wavelength of the spectral range. 
Let us assume that each subsystem of the ACS must make the same amount of change in the image quality of 


the OED and T,,(N) = const during observation. Then, the definition of acceptable MTF R, VPS, AFS, will be 


simplified (3): 
OT yV 
TP") = * Fac G = JLT,B, ). 


We expand functions T;(N), T-CN), Tg(NV), Tp CN) (3) into a series and obtain expressions that determine the TDE: 
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aaa 0,824 j1—7A°" Gn) _ 0,335 1-72") 
Ag™ < — +, Agi™ < = —1___ (5) 


, 
Vu Vu 


0.5 in[ 73 (va))-+ 
Ag S 3,018.71 —T3°"(v,), Aan” S _—— 
Here, v,, — Nyquist frequency. 
The limiting spatial frequency that the OED should resolve during the observation process is determined by the 


Johnson criterion: 


NyL 
va Tag? V3 =05vii (6) 
where v,,— angular (bar/rad) limiting spatial frequency; Ny — number of resolution elements (Johnson 


numbers); L — distance to OO; Ay, — critical size of OO. 

When decomposing, all members of the series are discarded, except the first ones, therefore, the tolerances obtained 
should be revised by gradient descent. 

The above reasoning allowed us to determine the accuracy of image stabilization of the onboard optoelectronic 
system (OOES) for tracking and stabilization modes. A methodology and a program for calculating MTF and tolerances 
were developed (Patent RU2021660340’). 

Let us denote requirements for the dynamic error of the ACS OED, taking into account the permissible MTF OED. 
Initial data for calculation: 

— diameter of the entrance pupil — D = 90 mm, 

— OS focal length — f*= 100 mm, 

— wavelength — i= 4 um, 

— distance to the object of observation — L = 10 km, 

— meteorological visibility — SM = 15 km, 

— flight altitude — H = 5 km, 

— number of MPD elements — 256X256, 

— frame rate — F = 400 Hz, 

— probability of correct detection — P 5, = 0.8, 

— false alarm probability — P,, = 10°, 

— signal-to-noise ratio of the amplifier — optimal filter — m = 1.105, 

— angular size of the object of observation — y; = 5-10 ° rad, 

— angular distance between two adjacent targets — y, = 3.34- 10 °rad. 

MPD parameters: 

— d=0. 0.25a — thickness of the photosensitive layer, 

— €, = 4:10 ° — charge transfer inefficiency, 

— Mc = 2 — number of samples per 1 element, 

— k, = 0.6 — fill factor, 

—~1=5.4-10 °s — reading time, 

— Tap = 5-10 '°s — time constant of the transducer. 

Figure 2 shows the MTF of linear, harmonic and vibrational shift of the image corresponding to the permissible and 
obtained DDP of the ACS OOED. 


' Program for calculating the modulation transfer function of an optoelectronic device based on UAV: RF patent 2021660340, 2021. Innopolis 
University. (In Russ.) 


Machine building and machine science 


153 


http://vestnik-donstu.ru 


Advanced Engineering Research 2022. V. 22, no. 2. P. 150-160. ISSN 2687-1653 


T 10 
0.8 
0.6 —— Tn. (Vv) 
== T rapa. (v) 
0.4 ° ——— Tansy. (V) 
0.2 
0. 


0 100 200 300 400 500 600 -Y, I/pan 


Fig. 2. Modulation transfer functions T°"(v) 


In accordance with (3)—(5) and the initial data for the calculation, we determined: v,, = 600 rad' of the OOED 
observation channel in the infrared region with a detection probability of P =0.8. For MTF OOED and TDE ACS: 
Ady = 3.5 ang. min, Aor = 3.1 ang. min, Aog = 3.0 ang. min. 

On the basis of the Lagrange equations of the second kind and the mixed Gilbert method, a mathematical model of 
the SO OED [4] of guidance and tracking was obtained in the form of three information channels driven by DBM 
electric motors along the gimbal axles. We presented it in a matrix form’: 

A(q)q + B(t,q)q + F(@.q) +W,q) + 0,9) = May — Myp. (7) 
Here, q = (a 8)", a, 8 —CO rotation angles in azimuth and elevation.. 


_ (Bi + B(B) —D(B) - 0 bag(t, a, B) 
aca) = ( —D(B) Gy ).20.0= (9 cerca) 0 ) 


. _ (—24BF (a, B) — B7E(B) _ (Walt, q) 
ae ( PF (B) } wed = (vce) 
r(q) ~ ~ _ Mapa _ Myp1sign(&) 
Ota) [Gp] OAR Ma = (42s) Mo = (mn) 
bog (t, a, B) = {(2A2(B) — Cz) cos a — 2E(B) sin A} Ox, (t) - 2F (B) Wy, (t) 
— {(242(B) — Cz) sin a + 2E(B) cos o}z,(t), 
Walt, q) = {-@x, (CF (a B) + @y, (t)(By + B(B)) — wz, (t)D(a, B)} + {E (a, B)NZ, (t) — E(a, B)w?, (t) + 
+ D(a, B)a(t)w, (t) + 24s (a, Box, (ez, (t) — F(a, Bay, (Ho2, (O}, 
we(t,q) = {we, (t)[E(B) cos (a) — Cz sin(a)] — &y, (t)D(B) + wz, (CE (B) sin(@) + C, cos(a)]} - 


i 
= 5 (Fea + cos(2a)) + D(B) sin(2a)]w%,,(t) — 2F (B)w5, (0) + [F(B)(1 — cos(2a)) — D(B) sin(2a)]wz, (t) + 
+2[24,(B) cos (a) — E(B) sin(a) Joy, (t)wy, (t) — 2[F(B) sin(2a) — D(B) cos(2a)]a,, (t)z, (t) — 

—2[E (B) cos(a) + 2A2(B) sin(a@)ay, (tmz, (t)}, 

—X¢, Sin(a) + Zc, cos(a) . 
Te (@) =m, 0 + 


—X¢, COS(M) — Zc, sin(a) 


—(X¢, cos(B) — yc, sin(B)) sin(a) + Zc, cos(a) " 


+Mz 0 
—(X¢, cos(B) — yc, sin(B)) cos(a) — Zc, sin(a) 


*Burdinov KA, et al. Matematicheskaya model' i sintez sistemy avtomaticheskogo upravleniya bortovogo optiko-ehlektronnogo pribora. In: Proc. XII 
All-Russian Congress on Fundamental Problems of Theoretical and Applied Mechanics, Ufa, 2019. P. 190-192. (In Russ.) 
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(-x¢, sin(B) — Yc, cos(B)) cos(a) " 
rg(q) = mz X¢, cos(B) — Yc, sin(B) 
(xc, sin(B) + Yc, cos(B)) sin(a) 
The next stage of development was the decomposition of doubly connected ACS [5] with nonlinear cross couplings 


in the CO, which is the studied OED. The linearized equations of motion with respect to the trajectory have the form: 
@ Xzy (t) = Oot € 0 + 27; B X,, (t), = Bot € 0 + 11/2, 
(Ayo + 44, )AG + (Bry + By )AG + Cy, A + a,2AB + by AB + Cy2AB = Ky Ary, 


Az, AG + by, AG + Az ,AB + bygAB + C22AB = Ky2Au2, (8) 
Au, = R,(p)Ag,, Ae, = Aaj, — Aa*; Aug = R2(p)Ae2, Ae2 = AB, — AB”. 


Having constructed the variable coefficients of system (8), we determine their values at ft = ¢, at “dangerous” points 
at which the parameters change significantly or change signs in the range of angles (B = 0 + 7/2, a =0+ 2n). Then, 
reducing system (8) (for t = t, ) to a system with constant parameters and using Cramer's rule, we express it in the form 
of transfer functions (TF) at points t = te: 

Aa = Wi; (p, ty Au, — Wy2(p, ty )Au2, (9) 
AB = W22(p, ty Aug — Wai (p, ty) Au. 

The obtained TF (9) taking into account (8) for each selected moment of time (t= ty ) can be presented in the form of 

a block diagram with direct cross couplings (CC, Fig. 3). From now on, for simplicity of writing, we denote TF as 


Weikp) = Wikpate )- 


Fig. 3. Structural diagram of doubly connected ACS 


We evaluate the CC ACS by their frequency response (FR). For this purpose, we write down the TF of the open 
system, taking into account the 2nd closed control channel when Af, = 0 (Fig. 4). 


Wi (p) 


—> 


Fig. 4. Structural diagram of ACS, open on one control channel 
Let us build hodograph wig ©) = Ry(jo)W,.41 (jo) and tubes around it with radii for each moment of time ft, : 
E1x(@) = Ri (®) Wai (@) Wx (@) 5x (@). (10) 
Taking into account (10) and the Nyquist criterion, it is possible to judge the stability of the ACS (1), (2) and the 


influence of the CC on its stability in the range of scanning angles (3). 


As a result, for each (t, -th), we determine the CC time: 
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Wy12(@) Wx21(o) oar 
8 xi (@) a 8x12 (0) 6,.21(@)W,; (o) = ae Waal W,;(@), jij a 1,2). (1 1) 


Let us consider the CC ACS that satisfy the decomposition conditions. 
According to (8), (11), from the CO mass geometry data, for 11 selected points in time, we determined (ty = 0; 0.89; 
1.35; 1.58; 1.71; 2.3; 2.97; 3.8; 4.5; 4.75; 6.28 s). We built a,(2), bj(t), c,(d). Logarithmic frequency response (LFR) 


5,(@) (kK = 1.11) was created in the Mathcad environment (Fig. 5). 


9 | — is MELEE 
-100 — k=] 
Tees —k=2 
9 -150 | (ee 
3 i 
B -200 | | pee 
< -250 —k=8 
—k=9 

-300 | 

10° 10° 10? 10* Frequency, rad/s 


Fig. 5. Logarithmic frequency response 6,(@) 

Figure 5 shows that the entire LFR set does not exceed —50 dB. From the analysis of LFR 5,(@), we have the results 
presented below. 

YoeQ, = (10 “+10*) rad/s. 

The following was determined: max 65,(@) = 0.0136 (-37.3dB)<8"=0.0335 (-29.5 dB), 
AL = M*°"/ (M" + 1) = 0.535 (5.43 dB) at Miop= May= 1.1; WM" = 1.15. 

Hence, a doubly connected control system (8) of the OED with the given parameters of the CO for the motion in 
question can be presented as two independent control channels in azimuth and elevation angle. 

Research Results. Figures 6 and 7 show the OOED control processes in stabilization and tracking modes. The scale 
of the misalignment graphs has been increased by 600 times. The graphs of the program control and the output values 


visually coincide. 


30 


20 


Angle, deg 
° 


0 l 2 3 4 5 6 7 


Time, s 
program, deg 
eSiemeee misalignment, deg x 600 


a output, deg 


Fig. 6. OOED stabilization processes in azimuth 
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Fig. 7. OOED tracking processes in azimuth 


We describe the harmonic oscillations of the carrier: 4; = A.= 2 (amplitudes), 7; = 7,= 1 s (oscillation periods), and 
A,=A,= 12 T, = T,=6.28 s. We indicate the load moments: M,,; = 2.2 Nm, M,,; =0.3 Nm. With such harmonic 
oscillations and load moments, the error of stabilization of the sighting axis does not exceed: 

— in azimuth — 1.2 ang. min, 

— in elevation angle — 1.6 ang. min. 

In tracking mode, with an additional control action introduced in the form of a constant speed of 12 desiy the error of 
tracking and stabilization does not exceed 1.7 ang. min in azimuth and 1.3 ang. min in elevation angle. 

To study the dynamics of the spatial control of the OOED in accordance with the ACS technique, the CSM of the 
digital automatic control system (DACS) of a specific OED was developed (Fig. 8) [1]. 


180/3.14°60, 


beta Regulator 


OUTPUT 


, [Ch 180/3.14] { 
Software 


device 


Au 
Lie}. . ges ‘ 


Carrier Object 


Fig. 8. Computer simulation model of a digital system for automatic 
control of an optoelectronic device 


The solution is implemented in the Matlab Simulink environment. Its elements: 

— CSM CO (Object) obtained according to (7); 

— CSM of drives (Drive); 

— CSM of PID-regulators (Regulator) taking into account non-linearities and DCD; 

— CSM of software guidance device; 

— Carrier CSM that implements the equations of carrier motion. 

Figures 6 and 7 show the results of modeling the DACS OED without taking into account the movement of the 
board. In this case: 

— guidance error does not exceed 18 ang. min in azimuth, 10 ang. min in elevation angle; 

— tracking error — 0.6 ang. min in azimuth, 0.6 ang. min in elevation angle; 

— stabilization error — 0.45 ang. min in azimuth, 0.6 ang. min in elevation angle. 

Figure 9 shows the control processes of the OED, taking into account its movement. The azimuth guidance error is 
2 ang. min, the elevation angle error is 5 ang. min. The error of tracking and stabilization is 0.8 ang. min in azimuth, 


1.3 ang. min in elevation angle. 
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Fig. 9. Simulation results of the system in two-channel guidance mode (degrees): 
a — system misalignment; b — transient schedules. In second case, schedules of program control (alpha target, beta target) and 
output value (alpha, beta) visually coincide 


Discussion and Conclusions. The results obtained are presented below. 

1. A technique for evaluating tolerances for the accuracy of image stabilization of the OOED under conditions of 
controlled and temperature influences was proposed. (Patent RU2021660340). Analytical estimates of the TDE image 
stabilization according to (4), (5) for tracking and stabilization modes based on the MTF OED were presented. The 
requirements for the ACS were obtained, providing a solution to the problem of detecting OO in the infrared region. 

2. The decomposition acceptance criterion was given, the method of decomposition of a doubly connected ACS 
OEP with the use of computer technologies was described. The requirements for the design parameters of the CO of a 
real device were defined. That allowed synthesizing the ACS OEP as two independent control channels. 

3. The CSM ACS OED was proposed to study the dynamics of isolated control channels and spatial control. In this 
case, the nonlinearity and non-stationarity of a specific CO were taken into account. 

4. A mechanical model of the OED was developed, the equations of the dynamics of spatial motion (7) of the CO 
were obtained. 

5. The parameters of the PID-controllers were synthesized, providing stable guidance and stabilization of the OED 
under all specified modes. Q-factor in speed: K,, = (3700 + 4000)s71, K,, = (3000 + 3300) s~?. Corrective links: 
T11 = Tez = 0.058, Ty21 = Tya27 = 5 X 10-5 s. This provided more accurate operation of the flight planning system 
and trajectory construction [6-8]. 

6. It was shown that with the help of CSM, it was possible to obtain the required time and error of guidance, 
stabilization, and tracking. In this case, the guidance speed was 0.6 s. The corresponding indicator of foreign 
analogues — | s°. 

7. In [9], a laser counteraction system was considered. Its technical characteristics allowed performing the task at 
W= 1000 W/steradian and laser power P=3.1 W. It was shown in[10] that the required radiation power was 
W = 2700 W/steradian at a laser power of P = 200 MW. 


> Directed Infrared Countermeasure (DIRCM). EMSOPEDIA. emsopedia.org. URL: https://www.emsopedia.org/entries/directed-infrared- 
countermeasure-dircm/ (accessed: 01.05.2022). 
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We calculated the required stabilization accuracy for each of the presented cases using formula 
2x180x60 


A arccos (1 = —) [ang. min] and made Table 1. 
Table 1 


Comparison of required stabilization accuracy 


Radiation power (W), | Laser power (P), | Required stabilization accuracy (A), 
W/sterad WwW ang.min 
Installation [11] 2,700 0.2 33.4 
Installation [10] 1,000 3.1 216.1 
Product 1.3 


Thus, the accuracy of stabilization of the proposed system was ten times higher than that required in papers [9, 10]. 
And this reduces the requirements for the convergence of the laser beam and the laser radiation power tenfold when 
developing the optical path of the product in question. 


List of abbreviations and conventions 
OES — optoelectronic system, 

OED — optoelectronic device, 

OO — object of observation, 

CSM — computer simulation model, 
ACS — automatic control system, 
MTF — modulation transfer function, 
TDE — tolerated dynamic error, 

ISS — image stabilization system, 
MPD — matrix photodetector. 
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